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Measuring the deformation of a ferrogel sphere in a homogeneous magnetic field 
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A sphere of a ferrogel is exposed to a homogeneous magnetic field. In accordance to theoretical 
predictions, it gets elongated along the field lines. The time-dependence of the elastic shear modu- 
lus causes the elongation to increase with time analogously to mechanic creep experiments, and the 
rapid excitation causes the sphere to vibrate. Both phenomena can be well described by a damped 
harmonic oscillator model. By comparing the elongation along the field with the contraction per- 
pendicular to it, we can calculate Poisson's ratio of the gel. The magnitude of the elongation is 
compared with the theoretical predictions for elastic spheres in homogeneous fields. 
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I. INTRODUCTION 

A sphere is the most perfect, most symmetric geomet- 
rical object. The Pythagoreans believed that the sun and 
the earth are perfect and thus spherical. Of course the 
rotation of the Earth is breaking the spherical symmetry. 
Whether this changes the shape of the Earth to a prolate 
or oblate rotational ellipsoid has been the subject of a 
long lasting quarrel between Newton and Cassini [l[ . The 
breaking of the spherical symmetry is found to be impor- 
tant in various fields of physics. The efficiency of nuclear 
fission depends on the shape of the nucleus Q. Partic- 
ularly convenient is a breaking of the symmetry by an 
external electric or magnetic field. A drop of a magnetic 
liquid, a colloidal dispersion of magnetic nanoparticles 
3 , elongates along the direction of the applied homoge- 
neous magnetic field 0, 0, @] . Very recently an amount 
of "quantum ferrofluid" , consisting of dipolar Cr- Atoms, 
has been found as well to elongate in a magnetic field 0] • 

All examples presented above have in common that 
they deal with fluid bodies, without any elasticity. In 
contrast smart nanocomposite polymer gels are kept 
together by their elastic polymer matrix. Especially fer- 
rogels @[l0] are a promising class for many applications, 
like soft actuators, magnetic valves, magnetoelastic mo- 
bile robots [ll|, [ljj, artificial muscles [l3|, or magnetic 
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controlled drug delivery [14|. All these applications are 



based on the magnetic deformation effect, a coupling be- 
tween the mechanic and magnetic degrees of freedom. 

The magnetic deformation effect is studied in its sim- 
plest form for a spherical gel body subjected to a homo- 
geneous magnetic field. An approximation for the result- 
ing deformation has been given in 1960 by Landau for 
the case of a dielectric elastic sphere 1151 and can read- 
ily be transferred to ferrogels [l6|, UtI Il8l| . The relative 
elongation e in this case was calculated as 
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where k = 1 /i5, G is the shear modulus and M is the 



FIG. 1: Scheme of the experimental setup. 



magnetization. Recently, the elongation has been recom- 
puted by Raikher [13] without constraining the shape to 
an ellipsoid. In this case, the elongation is expected to 
be 30 % larger. This effect has not yet been observed, 
a possible reason being that it is rather small for the 
large values of G, characteristic for most of the cova- 
lently cross linked polymer gels [jj 0|. In contrast, the 
elasticity of the new class of thermoreversiblc ferrogels 
[2fj| can reversibly be tuned via their temperature. In 
the present paper, we cast thermoreversiblc ferrogels in 
spherical samples and expose them to a uniform magnetic 
field, to test the above predictions. 



II. EXPERIMENT 
A. Setup 

Our experimental setup is sketched in Fig.[T] A ferrogel 
ball is immersed in a rectangular container, positioned on 
the common axis midway between two Helmholtz coils. 
For the empty Helmholtz pair of coils, the spatial homo- 
geneity is better than ±1%. This grade is valid within 
a cylinder of 1 cm in diameter and 14 cm in height ori- 
ented symmetrically around the center of the coils. The 
coils are powered by a current amplifier (fug electronic 
GmbH), which is controlled by a computer. The mag- 
netic system cannot follow the control signal immedi- 
ately. For a maximal jump height AB = 36 mT the 
field is reached after ts = 30 ms, as recorded by the Hall 
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FIG. 2: Image of the symmetrical (left) and distorted 
ball (right). The dotted line displays a fit with a circle (left) 
and an ellipse (right). 



probe (Group3-LPT-231) connected to a digital tcslamc- 
tcr (DTM 141). 

The temporal evolution of the ball shape is recorded 
using a highspeed camera, capable of taking 400 frames 
per second with a resolution of 768 x 768 pixels. Fig- 
ure[2]shows the original and distorted shape. The dotted 
curves stem from a fit of a circle to the edges in the image. 
The circle is located utilizing a normalized correlation 
technique [2l[, where the correlation between the gradi- 
ent of the image and a circle is maximized. The gradient 
is estimated by the Sobel operator [22||, and the circle is 
rasterized with anti-aliasing provided by a Gaussian fil- 
ter [23j]. This method is able to extract the radius and 
the coordinates of the center with subpixel resolution. 



B. Material 

The ball was prepared with a magnetic, thermo- 
plastic elastomer gel or in other words with a ther- 
moreversible ferrogel [2(j. Therefore, an ABA- type 
poly(styrenc-b-(cthylcnc-co-butylcnc)-b-styrenc) (SEBS) 
triblock copolymer (Kraton G-1650) was used as a gela- 
tor. As measured with size exclusion chromatography 
(SEC), Kraton G-1650 exhibits a molar mass of 99 000. 
The styrene content is 29 % (manufacturer information) . 
The gelator concentration is 3.5 w% per paraffin oil used. 
The ferrofluid was prepared in the classical way [2(j with 
24.5 w% magnetite per paraffin oil Finavestan A 50 B 
(Total). The shear viscosity at 20 °C of llmPas and 
the molar mass of 280 (manufacturer information) are 
lower than that of the earlier used paraffin oil Finaves- 
tan A 80 B [20}. Using A 80 B, the resulting ferrogels are 
too stiff to be remarkably deformed by the magnetic field. 
On the other hand, using the gelator Kraton G-1652 with 
a lower molar mass of 79 000 (SEC), the prepared ferro- 
gels are still softer, i. e. they are even less shape-retaining. 
Furthermore, they sweat out ferrofluid, slowly on stand- 
ing, faster in the magnetic field. 
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FIG. 3: Magnetization curve of the ball. 



C. Sample preparation 

Above its softening temperature at 45°C the mate- 
rial becomes a magnetic liquid. We produce the ferrogel 
sphere by casting the liquid into an aluminum mould at 
55°C. The mould consists of two parts with a spherical 
cavity, tightly screwed together, which is connected by 
a thin channel to a hopper mounted on top of the up- 
per part. To avoid any air bubbles in the final sphere, 
the mould is first filled with the liquefied material un- 
der vacuum (« ImPa). Then atmospheric pressure is 
applied, which compresses any low-pressure air bubbles. 
By repeating the process of varying the ambient pressure 
and under the influence of gravity, eventually all bubbles 
leave the mould via the hopper. Next we cool down the 
mould and separate the two parts of it. 

The ball is then immersed into water with a temper- 
ature of 24°C, where all experiments have been carried 
out. The water contains an amount of salt, so that the 
density of the liquid is only slightly less than that of the 
ball. This is necessary to reduce the influence of the grav- 
ity on the shape of the sample; the ball is so soft that, 
without buoyancy, it is deformed into an oblate ellipsoid 
under gravity. As an additional benefit we can easily 
measure the density of the gel by examining the density 
of the fluid. The density determined by this method is 
p= 1.085 ± 0.005 g/cm 3 . 



D. Magnetization 

Next, we measure the magnetization of the sphere 
for various fields, utilizing a fluxmetric magnetometer 
(Lakeshore, Model 480). The resulting magnetization 
curve is plotted in Fig. [3] The solid line displays an ap- 
proximation with the model presented in Ref. [241 ]. The 
sample is superparamagnetic with an initial susceptibil- 
ity xo = 0.81. 
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FIG. 4: The force response to a jump in the deformation. 
The solid line displays a fit to Eq. ([2} with the parameters 
Grheo = 1480 Pa, t = 127 s and /3 = 0.218. 
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FIG. 5: Elastic response of the ball, when a magnetic field is 
suddenly applied. The solid line represents a fit to Eq. ([3]). 



E. Rheology 

For the characterization of the elastic properties we 
measure the shear modulus G. with a rheometcr 
(MCR301, Anton Paar) in cone and plate geometry. The 
cone has a diameter of 50 mm and a base angle of 1°. Wo 
perform a stress relaxation experiment on the sample: 
from the equilibrium position we shear the sample by 
a deformation of 7 = 1 % and measure the stress r as 
a function of time. Figure [4] displays the results. The 
restoring force decays by 50 % during one second, which 
means that the shear modulus G = r/7 cannot be treated 
as constant. This decay can well be approximated by a 
stretched exponential (the solid line in Fig. 2]) 

G(t)=G lheo exp(-(t/t f). (2) 

The material therefore softens under load. To account for 
that time dependence, the magnetic experiment cannot 
be performed in a static manner. 



1.5 




10 20 30 40 50 60 70 80 90 
Magnetization M 2 (kA/m) 

(b) 



FIG. 6: a) The static elongation as a function of M 2 . The 
circles (triangles) denote the elongation (contraction) par- 
allel (perpendicular) to the applied field, respectively. The 
solid (dashed) line is the best linear fit. b) The frequencies of 
the initial vibrations obtained from Eq. ©. 

III. RESULTS AND DISCUSSION 

We performed time-resolved measurements of the rel- 
ative elongation e of the ball after applying a magnetic 
induction B in a jumplike fashion for ten different val- 
ues of B. Figure [5] shows e(t) for the case B = 36 mT. 
The elongation e = (d — d )/d measures the scaled dif- 
ference of the diameters in the direction of the field (d) 
and without a field (do). Due to the sudden increase of 
the magnetic field and the inertia of the ball, the latter 
performs uniaxial damped vibrations with a frequency 
/ = 22.4 Hz. When the oscillations cease, the elongation 
continues to grow due to the softening under load. Be- 
cause the timescalc of the softening is much larger than 
that of the oscillation, the experiment may be approxi- 
mated by a harmonic oscillator pulled by a constant force, 
where the spring constant relaxes according to Eq. ^ 

e + Se + ojI cxp (-(*/*„)") e = F. (3) 

Here S is the damping constant, F the pulling force (i.e. 
related to the magnetic induction) and too the natural fre- 
quency. A fit with that equation is displayed in Fig. [5] as 
the solid line. The ratio F/ljq corresponds to the relative 
elongation in the equilibrium state if the spring constant 
would not change. Therefore, this ratio compares to the 
elongation predicted by the static theories. 

The dependence of the elongation on the magnetiza- 
tion is shown in Fig. [S]a. For ten different values of the 
magnetic induction we have recorded and evaluated the 
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elongation e of the ball. Figure [6]a presents the outcome 
for the elongation parallel (e z ) and perpendicular (e x ) 
to the magnetic field. The data have been plotted ver- 
sus M 2 . In agreement with Eq. (JT|) we find a linear re- 
lationship Si = c t M 2 with c z = 13.4 • 10 _5 (||) 2 and 
c x = —6.1 • 1CP 5 (j^-) 2 . Figure 0b shows the oscillation 
frequency / = uiq/(2it) versus M 2 . 

Next we compare the ratio of the slopes c z /c x = k z /k x 
with the theoretical predictions. For a deformation re- 
stricted to an ellipsoidal shape and the assumption of a 
uniform strain field, the expressions given in [pa . Il7| can 
be rewritten in terms of Poisson's ratio a 
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For an incompressible gel (a — 1/2) this leads to k z = jg. 
For arbitrary a one obtains for the ratio 
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We exploit Eq. to determine a, and by substituting 
a in Eqs. (|4bll[) one finally arrives at G. This yields 



(7 = 0.48 + 0.01 G = 0.65kPa. 



(6) 



For the more general case of a non-uniform strain field 
and a shape not restricted to an ellipsoid [Tt} , one obtains 
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which yields 
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and finally gives 



(7 = 0.47 + 0.01 G = 0.87kPa. 



(7a) 
(7b) 

(8) 
(9) 



Both approaches yield a Poisson's ratio a very close 
to the limit of incompressibility a = which is char- 
acteristic for rubberlike materials [25|. As reported be- 
fore , the values derived for the shear modulus G differ 
by w 30%. But the value obtained from the rheometcr, 
G r hoo = 1.48 kP a, exceeds the old and new predictions 
by a factor of 2.3 and 1.7, respectively. So none of them 
is corroborated by the experiment. 

The frequency of the vibrations after the sudden in- 
crease of the magnetic field offers another possibility to 
measure the shear modulus. For an incompressible clas- 
tic sphere that performs spheroidal vibrations, where the 
sphere gets alternately deformed into a prolate and oblate 
ellipsoid of revolution, the frequency is given [26[ by 
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FIG. 7: The shear modulus measured by different methods: 
Commercial rheometer (dashed), elongation theory [ItJ (short 
dashed), elongation theory [l5T| (dash-dotted), vibrations (solid 
circles) . 
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FIG. 8: The shear deformation as a function of time, 
applying a constant deformation of 1.2% for 800 s. 
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Since the radius of the sphere r and the density p are 
known, we can compute the shear modulus from the mea- 
sured vibration frequency /. Figure [7] shows a compar- 
ison of the value of the shear modulus obtained from / 
with the values from the elongation and the mechanical 
measurement. The average shear modulus determined by 
this method is G v ib = 0.1 kPa, which differs by a factor 
of 15 from G r hco- This large deviation may arise, because 
the model for the vibrations does not take into account 
the surrounding water. This needs to oscillate together 
with the sphere, leading to an increased effective mass of 
the oscillator, and thus a reduced frequency. 

Now we reconsider the deviation between the shear 
modulus G as determined from the elongations and G r h e o 
determined via rheology. They differ by factor of about 
two. This difference cannot be explained solely by the 
inaccuracy of the measurement of G r heo, which is typi- 
cally around 10 % for a commercial rheometer. Rather, 
the deviation is likely to stem from the fact that the ide- 
alistic models [TH, [TtJ consider only Hookean elasticity, 
which is insufficient for our material. In fact, the defor- 
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mation of a gel put under load can be made up of three 
contributions, namely instantaneous elastic deformation, 
retarded anelastic deformation and viscous flow [13] ■ Our 
material clearly shows anelasticity, as demonstrated in 
Fig. |4j Moreover, our sample shows the phenomenon of 
viscous flow, as illustrated in Fig. [H In this measure- 
ment we first apply a constant strain for 800 s, and then 
record the relaxation of the strain for zero stress. The 
remaining deformation at t = 7000 s is about 15 % of the 
strain initially applied. This value can be regarded as an 
upper bound for the plastic contribution to the deforma- 
tion (viscous flow). While the existence of anelasticity 
and viscous flow gives no straight-forward explanation 
of the 70 % deviation between experiment and theory, it 
is obvious that the full viscoelastic behaviour of the gel 
should be included in the computations from the very 
beginning. 

IV. CONCLUSION 

We have measured the deformation of a ferrogel sphere 
in response to a uniform magnetic field by direct optical 
means. Wc compare the results for the first time with 
the models in Refs. fl5l. [l7|. From the ratio of the elon- 



gation parallel and perpendicular to the field, we calcu- 
late Poisson's ratio, which is close to the value a = 0.5 
expected for an incompressible material. More impor- 
tantly, the absolute value of the elongation is 70 % larger 
than the one calculated from the models. This is pre- 
sumable caused by the neglection of the anelasticity and 
viscous flow of the ferrogel. So further theoretical inves- 
tigation, including the full viscoelastic properties of the 
gel, is needed to explain the experiment quantitatively. 
With such a model one would also be able to compute 
the dynamic response of the gel under a sudden change 
of the external magnetic field. This is not only of funda- 
mental interest, but also important for possible technical 
applications of these smart materials. 
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